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Abstract 

Time-resolved magneto-optical Kerr effect (MOKE) and Brillouin light scattering (BLS) spec- 
troscopy are important techniques for the investigation of magnetization dynamics. Within this 
article, we calculate analytically the MOKE and BLS signals from prototypical spin-wave modes 
in the ferromagnetic layer. The reliability of the analytical expressions is confirmed by optically 
exact numerical calculations. Finally, we discuss the dependence of the MOKE and BLS signals 
on the ferromagnetic layer thickness. 
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I. INTRODUCTION 



Optical techniques based on the magneto-optical Kerr effect (MOKE), such as time- 
resolved MOKE (TR-MOKE) or Brillouin light scattering (BLS) spectroscopy are routinely 
used for the investigation of magnetization dynamics. Whereas TR-MOKE provides informa- 
tion on the magnetization dynamics within the time domain, BLS informs on magnetization 
dynamics within the frequency domain. 

BLS measures light intensity scattered by a spin wave. In a phenomenological (wave) 
picture a spin wave is a periodical displacement of the magnetization with respect to the 
saturated state, oscillating at a spin- wave frequency Wgw, and propagating at a given spin- 
wave A;- vector. The spin waves result from the coupling between the spins, dominated by 
exchange and dipolar interactions [IHS]. Any periodical variation of the optical properties 
(namely, periodic variation of the permittivity tensor e{t, f)) works as an effective oscillating 
and propagating optical grating |lH7j. Hence, due to the spatial periodicity of the effective 
grating, the reflected light is scattered and changes its propagation direction. Furthermore, 
the oscillation of the effective grating changes the light frequency of the scattered light, a 
quantity detected in a BLS spectrometer |H]. In a (pseudo-)particle picture, the scattering of 
the light by the spin wave can be interpreted as inelastic scattering of photons on magnons 
[9], where the photon is gaining (loosing) its energy as it absorbes (creates) a magnon, 
respectively. Therefore, in contrast to TR-MOKE, BLS can also detect non-coherent spin 
waves such as thermal spin waves. Note, that the BLS technique usually provides a larger 
experimental sensitivity when compared to the TR-MOKE technique. For example, it has 
been demonstrated to detect thermal spin waves on a Co monolayer jTO] . 

TR-MOKE investigations are based on repeated excitations of magnetization dynamics, 
usually by magnetic field pulses or by intense light pulses pTHT3] . The resulting excitations 
are then stroboscopically detected making use of the MOKE. Therefore, TR-MOKE pro- 
vides an insight into the magnetization dynamics within the time domain. Using Fourier 
transformation, the TR-MOKE signal can be easily transformed to the frequency domain 
|14j . Hence, in case of externally excited systems, TR-MOKE and BLS investigations are 
complementary. This was nicely demonstrated by K. Perzlmaier et al investigating confined 
spin- wave modes in a square permalloy element [15]. 

The numerical models to calculate the BLS light intensity were elaborated by J.R. 
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Cochran et al IH [5], followed by L. Giovannini et al [6]. Later, a simple relation between 
the complex Kerr angle and the BLS intensity was expressed by M. Buchmeier et al [7]. 
However, all those treatments of the BLS intensity are numerical ones, and up to now, there 
exists no analytical expression of the MOKE and BLS signals originating from spin-wave 
modes. 

Within this article, we present the analytical dependence of the TR-MOKE and BLS 
signals for several types of spin-wave modes. Such calculations can serve either for separation 
of MOKE and BLS signals from a single ferromagnetic (FM) material in a stack of FM layers 
or for the quantitative determination of the energy carried by each spin-wave mode. 
The article is organized as follows: In Sec. [IT] we establish a relation between the BLS 



intensity and strength of the MOKE effect. In Sec. Ill the analytical expression of the 
MOKE depth sensitivity function is developed, and we discuss its validity and properties. 



Section IV provides the analytical expressions of the MOKE and BLS signals. Finally 
Section |V] compares analytical expressions with optically exact numerical models, and we 
discuss in detail the FM-layer thickness dependence of the MOKE and BLS signals. 



II. RELATION BETWEEN MOKE AND BLS SIGNALS 

As the light travels through the FM layer, the light intensity is attenuated and the phase 
of the light is delayed. Therefore, sublayers of the FM material situated at different depths of 
the FM layer provide a different MOKE response to a given magnetization state. Due to the 
fact that the MOKE is linear in magnetization, the result can be written as a superposition 
of the single contributions coming from the different depths of the FM layer [TTf 1201422] : 




where $s = Vps/rss and $p = — Tsp/rpp are the s- and p- complex Kerr angles, arising 
from the FM film when the incident light is s and p polarized, respectively. The terms 
Txy, x,y = {s, p} stand for components of the reflection matrix. Ls/p{z) and Ps/p{z) are the 
complex MOKE depth sensitivity functions related to longitudinal and polar magnetization, 
respectively, d being the FM layer thickness. mL{z,t) and mp{z,t) are the depth profiles 
of the magnetizations in the FM film having longitudinal (i.e. in-plane and parallel to the 
plane of the light incidence) and polar (i.e. normal) directions. In the case of spin waves. 
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those magnetizations correspond to profiles of tlie dynamic magnetization (i.e. spin-wave 
amplitudes), precessing at the frequency Ug^, 

mp{z, t) = mp{z) sinuswt (2) 

rriLiz, t) = mL{z) cosUswt- (3) 

According to Refs. [3 the BLS intensity of the backscattered light can be 

expressed in a rather similar way as the complex Kerr angle 

2 

(4) 



.d 2 







Tss/pp [-Ls/p{z)mL{z) + Ps/p{z)mp{z)] dz . 

Within the depth sensitivity, the main difference between the complex Kerr angle $ and the 
BLS intensity jg given by the fact that whereas MOKE is a linear combination of 

Lg/p, Pg/p and the magnetization profiles niL^z), mp{z), respectively, the BLS intensity has 
a quadratic form. 

From comparing the equations expressing the complex Kerr angle (Eq. [T| and the BLS 
intensity (Eq. |4]), their close similarity is apparent. With exception of the sign of the longi- 
tudinal contribution, the BLS intensity is basically a quadratic form of MOKE. Therefore, 
the BLS intensity can be expressed as being proportional to the square of the off-diagonal 
reflection coefficients r^p/ps 

^(.0±-.w) ^ |^^^^^^|2 ^ ^ (5) 

where we must reverse the sign of the longitudinal contribution when expressing $s/p- For 
example, this can be achieved either by reversing the sign of mi in the calculations. The 
complex Kerr angle (i.e. the magnitude of the MOKE) is denoted by $s/p, neglecting its 
time dependence (i.e. omitting the term sin(ci;swit + 0s/p) in Eq- ([l|)- 

III. ANALYTICAL EXPRESSION OF THE MOKE DEPTH SENSITIVITY 
FUNCTION 

In the case of an optically thick FM layer (i.e. the FM film thickness d is larger than the 
MOKE probing depth Amoke = A/(47rIm(A^^)), A being the vacuum light wavelength and 
Nz defined just below), Ls/-p{z) and Ps/p{z) can be analytically expressed as |2T] 

P./p{z) = P3/p(0) exp[-4mN,z/X\ (6) 
L,/p{z) = L3/p(0) exp[-AmN,z/\], (7) 
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where we define 7s/p to be the ratio of the LMOKE and PMOKE response at the upper 
interface (i.e. z = 0) of the FM layer, Ls/p(0) = 7s/p-Ps/p(0), is the normalized wave-vector 
of hght in polar (i.e. normal) direction, = \/ (iV(^™))2 — (A^('^i'') sinf/))^, where iV*^^™^ and 
are the refractive indices of the FM layer and air, respectively, and is the angle of 
light incidence with respect to the sample normal, respectively. In general, metals provide 
a relatively large value of the optical permeability eo = A^^, being the refractive index. 
Therefore, the angular dependence of A^^ can be neglected and hence A^^ ~ A^. For example, 
for Ni at A = 810 nm, — — 13.24 + llSSli and hence A^^ is reduced only by 1% when 
going from = to 90° [17j. 

An example of L^j^^z) and P^i^{z) is shown in Fig. [TJ^a) for the multilayer structure 
air/Cu(2nm)/Ni(60nm)/Cu(5nm)/Si. The calculations were done using a 4 x 4 matrix 
formalism [2l], for a light wavelength of A = 810 nm and an incidence angle of = 25°. The 
diagonal permittivity of the used materials are ^o^"^ = — 26.37+2.61i, eg'^^'' = —13.24+22.072, 
4^'^ = 13.58 + 0.04z, the off-diagonal permittivity of Ni being ef^^ = 0.217 - 0.091i [23]. 
As the difference between and Lp (Pg and Pp) is only the starting amplitude and phase 
of Ls/p(0) and Ps/p(0), the calculations in Fig. [l](a) are presented only for Ls and Ps. In 
this particular case, the polar MOKE (PMOKE) amplitude is about 12 x larger than the 
longitudinal MOKE (LMOKE) amplitude. It can be considered a general rule that the 
PMOKE is stronger than the LMOKE. Additionally, in our example, also the small incidence 
angle of (f = 25° contributes to the small value of the LMOKE (LMOKE vanishes at (f = 0). 
But even for an incidence angle of about 60-70°, when the LMOKE reaches its maximum, 
its amplitude would increase only by a factor of 2, still much smaller than the PMOKE. 

Figure [l](b) shows an agreement between the analytical expressions of the MOKE depth 
sensitivity functions Ls/p{z), Ps/piz) as given by Eqs. (|6]-[7]) and as determined from optically 
exact 4x4 matrix calculations for various Ni thicknesses. The starting point, Ls/p(0), Ps/p(0) 
is determined by the optical 4x4 calculations in both cases. It is demonstrated that there 
is nearly a perfect agreement for large Ni thicknesses 40 and 60 nm, whereas there is a 
disagreement for Ni thicknesses below 30 nm. It is because the analytical expressions of 
Ls/p{z) and Ps/p{z) are valid when the thickness of the FM layer d is larger than the MOKE 
probing length, Amoke = A/(47rIm(A^2)) (Amoke = 14.5 nm for Ni at A = 810 nm). On 
the other hand, for very small thicknesses of the FM layer (when d <^ Amoke? i-e. below 
d = 3nm), the FM layer can be neglected from the optical point of view. Hence, Ps/p, L^/p 
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can be considered constant so that the expressions of L^/p, Ps/p are also vahd for a very 
small thicknesses of the FM layer (so-called ultrathin FM layer approximation |26j). 

Therefore, the MOKE depth sensitivity functions Ps/p, L^/p are described well by the 
analytical expressions (Eqs. (|6]-[7])) with exception of the FM thickness range between about 
3 nm to 30 nm. However, even in this range, the analytical expression describes basic features 
of the depth sensitivity functions. Namely a reduction of |-Ps/p|) l-^s/pl and an increase of 
the phase arg(Ps/p), arg(Ls/p) with increasing depth inside FM film. Hence, deviation of the 
simple analytical calculations from an optically exact calculations is not large even in case 
of this thickness interval, as shown later. 

We finally note that Ls/p{z), Ps/p (-2) are nearly independent of the incidence angle as 
their angular dependence is governed solely by A^^, whose angular dependence is very weak 
in the general case of metals. 



IV. MOKE AND BLS FROM DAMON-ESHBACH AND PERPENDICULAR 
STANDING SPIN- WAVE MODES 

To obtain the MOKE or BLS response of a given spin-wave mode, the profile of the 
dynamic magnetizations mi{z), mp{z) through the FM film must be determined first. Those 
calculations are usually based on phenomenological models of the magnetization inside the 
FM layer P [3, ETJ EH]. Then, the complex Kerr angle $ or the BLS intensity /("^oi"^-) 
coming from a given spin- wave mode can be expressed using Eq. ([T]) and Q, respectively. 

In general, the profile of a spin-wave mode must be calculated numerically. However, 
prototypical spin-waves modes have rather simple analytical expressions of their amplitude 
profiles. Here we work out the magneto-optical response of three spin-wave modes, the 
Damon-Eshbach (DE) mode (including the homogeneous ferromagnetic resonance (FMR) 
mode) bounded to the upper (Fig. [2|^a)) and lower (Fig. |2|^b)) interface, and the perpendic- 
ular standing spin- wave (PSSW) mode (Fig.|2{c)). 

The DE mode occurs when M lies in-plane and the spin wave propagates in a direction 
perpendicular to M. This mode can be bounded either to the upper (Fig. [2|a)) or lower 
(Fig. [2|^b)) interface, depending on the mutual direction of the saturation magnetization and 
/c- vector propagation [29]. In case of DE mode bounded to upper interface (Fig. [2]^a)), the 
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polar and longitudinal profiles of the dynamic magnetizations are [29] 

m^^^^ {z, t) = m-p exp(-/c2^sw-2) sin oj^^t (8) 
j^(DEi) _ ^{o)g exp(— /c2_sw-2) cos uj^^t (9) 

where e = m)^^ /m~p^ describes the ellipticity of the precessing magnetization, assumed to be 
constant over the whole FM layer thickness, fc^^sw is the normal direction of the spin-wave 
wavevector. In the case of an homogeneous FMR mode, kz^sw = 0. Substituting the dynamic 
magnetization profiles (Eqs.([8||9|) and the MOKE depth sensitivity functions (Eqs. (|6|j7]) 
into the expression of MOKE effect (Eq.|l|), we get 

*s/p id,t) = myP,/p{0)Jl + 72/pe2 — — — sm{ujs^t + 0s/p) (10) 

* '^z,sw ~r lot I (X 

where d is the thickness of the FM layer and a = ArcNzd/X. As the detected MOKE signal is 
a mixture of both LMOKE and PMOKE, it results in a phase shift (ps/p between the phase 
of the spin-wave mode and the detected MOKE signal, tan0s/p = 7s/pe. Moreover, the term 



1 + Ts/pC^ originates from Pythagorean sum of PMOKE and LMOKE. 
In the case of a DE mode bounded to the lower FM interface (Fig. [2|b)), the profiles of 
dynamic magnetizations are analogous to Eqs. (|8]l9]) 



m 
m 



^p^^\z,t) = mP exp{-kz,sw{d - z)) sinUsy^T (11) 
^°^^^(z,t) = mPeexp{-kz,sw{d - z)) coscUs^r (12) 



leading to the MOKE effect 

^(DE2), , ,^ (o)d ~2 ^ exp(-m) - exp(-A;^,swC^) . . - n n Q^ 

The last type of the spin- wave mode to be discussed here is the PSSW modes (Fig. [2|^c)), 
described approximately as a cosine function with its maxima pinned at the FM interfaces. 
Then, the amplitudes of dynamic magnetizations are 

'm^p^^^\z,t) = m!p cos{m7rz/d) sincUs^r (14) 
m^^^^'^\z,t) = mp^ecos{m7r z/d) cos UswT (15) 

where integer m denotes the mode number of a given PSSW mode. Substituting those 
magnetization profiles (Eqs. (14 15)) to Eq. ([T]), leads to the MOKE effect 



K/p id,t) = myP3/p(0) Jl + 7s/pe' oo _ o sm(a;swt + (t^./p)- (16) 



Please recall, that the analytical expression of the BLS intensity of different 

modes can be easily obtained from the complex Kerr angle <l>s/p using Eq. (5), — 



s/p 



-^0 1 '^s/p'^ss/pp I 



, where $s/p is given by Eqs. (10), (13), (16), and where the sign of LMOKE 



must be reversed, e.g. by reversing the sign of 7s/p or rriL. 



V. MOKE AND BLS SIGNAL FROM SPIN- WAVES IN A Ni FILM 

Figure |3] shows the MOKE effect from various spin- wave modes as a function of the Ni 
film thickness in air/Cu(2nm)/Ni((i)/Cu(5nm)/Si structure. The investigated modes are 
the DE mode, bounded to both the upper (DEI) and the lower (DE2) interface, expressed 
for kz^sw = lO'^m^^. Furthermore, the investigated modes are the FMR mode and several 
PSSW modes with different mode numbers. Within those calculations, for demonstration 
purposes, we keep the spin-wave ellipticity e constant, although obviously the spin-wave 
ellipticity is different for different spin-wave modes. 

(i) First, the calculations demonstrate a nice agreement between the exact optical cal- 



culations based on 4 x 4 matrix formalism and analytical formulae (Eqs. (10), (13), (16)). 
The largest disagreement is for Ni thicknesses in the range of about 3 - 30 nm, as for this 
range the analytical expressions of Ls/p{z), Ps/p{z), as given by Eqs. (|6])-([7]), are not exact, 



as already discussed in SecTII 



(ii) As the Ni thickness is increasing, the MOKE signals increase and then saturate. In 
Fig. |3| the saturation is clearly visible only for the FMR, DEI and PSSWl modes. The 
saturation roughly appears when the depth of the unique sign of the dynamic magnetization 
corresponds to the length 2Amoke- In the case of the FMR or DEI mode, the MOKE signal 
saturates roughly at 2Amoke ~ d- Within our Ni example, 2Amoke ~ 30 nm and hence 
the MOKE signal saturates roughly at d = 30 nm for the FMR/DEl mode, in agreement 
with Fig. [3] For PSSW modes, a unique sign of the dynamic magnetization holds for 1/4 
of the perpendicular spin wavelength, which is Asw,± ~ 2d/m (m is integer denoting the 
PSSW mode order). Consequently, the PSSW modes saturate when (l/4)Asw,± = d/{2m) ^ 
2Amoke- Hence the estimated PSSW mode saturation appears aX d = 4mAMOKE, which 
corresponds to 60, 120, 180 nm, etc. for PSSW modes with m = 1, 2, 3, etc., respectively. 
This is also in qualitative agreement with the calculations shown in Fig. [3j 

(iii) For a large thickness of the FM layer {d — t- 00), the MOKE signal 
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from DEl/FMR and PSSW modes saturate to the same value, being ^s/p{d,t) = 
mPPs/p{0)^l + jl^^e^iX/AniN,) sin(a;s„t + 0s/p), as follow from Eqs. (10), (16). The un- 
derlying physics is that the upper part of the FM layer, which is probed by light, has nearly 
constant amplitude of the dynamic magnetization. Hence, obviously, in case of PSSW modes, 
such saturation appears for larger thicknesses of the FM layer, as compared to DEl/FMR 
modes. 

(iv) For smaller Ni thicknesses, MOKE signals for PSSW modes are strongly reduced with 
increasing number of the PSSW modes. The reason is analogous to the discussion in point 
(ii): when the depth profile of the PSSW mode oscillates on a short distance as compared to 
the probing depth of MOKE, Amoke, then MOKE contributions from various depths of the 
Ni film cancel each other. Quantitatively, when 1/4 of the perpendicular spin- wavelength, 
(l/4)Asw,± = d/{2m), is shorter than half of the MOKE probing length, Amoke/2, then the 
provided MOKE signal is significantly reduced. Substituting Amoke = 15 nm, the reduction 
of MOKE signal is for Ni thicknesses below 15, 30, 45 nm for PSSWl, PSSW2, PSSW3, 
respectively, in agreement with Fig. [3] 

Figure |4] demonstrates MOKE and BLS signals originating from the spin waves inside the 
Ni film, where the spin-wave profiles and frequencies are calculated using phenomenological 
models [271 [2H]- The used magnetic properties of Ni are: saturation magnetization HqMs = 
659 mT (i.e. Ms = 520kA/m), exchange constant A = 4.7pJ/m (i.e. exchange stiffness 
D = 2.46meVnm^), Lande (7-factor g = 2.1, out-of-plane anisotropy = OkJ/m^ [HO] . 
The in-plane spin-wave wavevector is gy = 0, and the external magnetic field fioH = 50 mT 
is applied parallel to the plane of incidence. The optical parameters of the Ni film as well 



as its optical surrounding are described above in Sec. |III[ The dependence of the MOKE 
signals on the Ni thickness and the related spin- wave frequencies are presented in Fig. |4] for 
the FMR mode as well as for several PSSW modes. The MOKE signals are expressed as 
|$|, Kerr rotation 6 = Re($) and Kerr ellipticity e = Im($). The amplitude of the in-plane 
dynamic magnetization of the FMR mode is chosen to be one, m^™^^ = 1. The amplitudes 
of the dynamic magnetizations of the other spin wave modes are normalized in a way that 
their energies per unit area are equal P^, 7] . 

All expressions of the MOKE signals (|$|, 6, e) provide very similar dependencies. As the 
external field H was applied parallel to the plane of incidence, the longitudinal contribution 
to the MOKE effect is zero, as there is no dynamical magnetization in the longitudinal 
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direction. However, as discussed in Sec. the longitudinal depth sensitivity function 
is here about 12-times smaller than the polar one. As the normal spin-wave amplitudes 
are about 1-2-times smaller than the in-plane amplitudes, the polar contribution would be 
dominant even in the case of H perpendicular to the plane of incidence, where the LMOKE 
contributes to the outgoing MOKE signal. 

The MOKE signals in Fig. |4] are compared to the normalized BLS intensity (dashed 
magenta line). As the BLS signal is basically square of the MOKE signal (Eq. ([s])), the 
behaviour of scaled BLS and MOKE signals are very similar, while the BLS signal is reduced 
more significantly for smaller Ni thicknesses or for higher orders of the PSSW modes. 



VI. CONCLUSION 



In conclusion, we have analytically expressed the MOKE and BLS signals originating 
from prototypical spin-wave modes, namely the Damon-Eshbach (DE) and Perpendicular 
Standing Spin Wave (PSSW) modes. The calculations are based on the additivity of the 
MOKE effect, on analytical expression of the depth sensitivity function of the MOKE signal, 
as well as on a straightforward relation between the MOKE and BLS signals. As a showcase, 
we have expressed the MOKE and BLS signal in a Ni film. It is shown that analytical 
calculations describe well the physical behavior, as follows from the comparison with exact 
magneto-optical calculations. Furthermore, we have demonstrated that the dependence 
on the FM-layer thickness of both the MOKE and BLS signals is very similar. Namely, 
with increasing FM layer thickness, the MOKE and BLS signals saturate, where for PSSW 
modes with higher mode number, the saturation is provided for larger FM layer thicknesses. 
Furthermore, the MOKE and BLS signals reduce significantly as the PSSW mode number 
is increasing. 
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FIG. 1: (color) (a) Depth dependence of the longitudinal and the polar depth sensitivity function 
Lg, Ps inside the 60nm-thick Ni film, (b) Depth sensitivity function P, for Ni-thicknesses 10, 20, 40 
and 60 nm, presented in the complex plane calculated by the 4x4 matrix formalism (dashed black 
line) and using the analytical expression Eq. ^ (solid red line). Light wavelength and incidence 
angle are 810 nm and 25°, respectively. Green arrows point the start point of the depth sensitivity 
function, Ps(0), i.e. its value at z = 0. 
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FIG. 2: (color) (solid red line) The sketches of the treated spin-wave modes (a) DE mode bounded 
to the upper interface (called DEI) (b) DE mode bounded to the lower interface (called DE2) (c) 
PSSW mode, (blue dashed line) Sketch of the MOKE depth sensitivity function. 
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FIG. 3: (color) MOKE signals |$| originating from the different types of spin waves as a function 
of the Ni thickness, assuming the spin- wave ellipticity e = 2. (t, A) denotes for the DE modes 
bounded to upper and lower EM interface, respectively, with fc^^sw = lO'^m^^. (★) denotes for the 
EMR mode (i.e. as the DE modes with /c^^sw = 0). (■) denotes for the PSSW modes. Incidence 
angle is 25°, light wavelength is 810 nm. Eor values of optical and magneto-optical parameters, see 
text. Symbols are the optically exact 4x4 matrix formalism, solid lines the analytical expressions 



(Eqs. (10 13, 16)). 
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FIG. 4: (color) Calculated MOKE signals (height of the surface plot: |<I>|, solid red line: Kerr 
rotation 9 = Re(<I>), solid green line: Kerr ellipticity e = Im(<I>)) and scaled BLS signal (dashed 
magenta line) as a function of the Ni thickness and the related spin-wave frequencies. For values of 
optical and magnetic parameters, see text. The top plane of the 3D graph presents the spin- wave 
dispersion of the Ni film (solid lines). The profiles of the spin- wave amplitudes through a 60- 
nm-thick Ni film are sketched above the 3D graph for four examples of the spin- wave modes, one 
FMR and three PSSW modes. Within those profiles, the left part is showing the trajectory of 
the end-point of the magnetization vector. The right part shows the amplitudes of the dynamic 
in-plane (solid line) and normal (dashed line) magnetizations. 
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